JOURNAL OF RESEARCH of the National Bureau of Standards- A. Physics and Chemistry 
Vol. 76A, No. 4, July-August 1972 

Anatomy and Thermal History 
of Laser Self-Focusing Damage Tracks in Glass 

Wolfgang Haller and Joseph H. Simmons 

Institute for Materials Research, National Bureau of Standards, Washington, D.C. 20234 

March 16, 1972 

Laser induced self-focusing damage tracks in glass are analyzed. The damaged regions consist of 
long cylinders of altered refractive index, containing occasional bubbles and partly surrounded by 
arrays of small cracks. An analysis of the refractive index change in the subfracture damage tracks and 
the presence of bubbles leads to estimates of a minimum temperature reached by heating of the glass in 
the tracks. A calculation of minimum stresses developed during laser-beam heating gives a good account 
of the orientation of the observed cracks. 
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1. Introduction 

Self-focusing damages caused by laser beams in 
glasses range from gross fractures ("starbursts" and 
planar cracks) to lines made up of very small bubbles 
("bubble tracks") [l]. 1 Less conspicuous permanent 
changes in the beam regions are the "filaments," 
variously described as bright or dark lines of 0.5 to 
2.0 /xm diameter [2]. These filaments appear as regions 
of altered refractive index and the discrepancies in 
their reported dimensions stem from difficulties inher- 
ent in the interpretation of microscopic refraction 
images. 

Figure 1 shows a typical damage track produced in a 
piece of glass. In picture B the objective lens focuses 
upon the track. In pictures A and C the objective lens 
focuses 9 /xm respectively behind and before the 
damage track. In the upper and lower parts of the 
picture, bubble tracks can be seen. These are separated 
by a bubble-free region in the center of the picture. A 
fracture in the glass appears near the lower bubble 
track. Pictures A and C show the refraction effect of a 
transparent, cylindrical region in the glass, which 
obviously has a different refractive index from the 
surrounding glass matrix. This region extends over the 
full length of the track and is the object of our present 
investigation. 

2. Production of Damage 

A 15 ns pulse from a nominal 10 MW ruby laser was 
focused into a polished cube of glass by means of a 
lens with a focal length of 9 cm. The glass had a com- 



position of approximately 61 percent Si0 2 ; 20 percent 
AI2O3; 14 percent Na 2 0; 4 percent Li 2 and 1 percent 
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1 Figures in brackets indicate the literature references at the end of this paper. 



0.1 mm 

FIGURE 1. Micrographs of a damage track illuminated in transmitted 

light. 

The objective lens is focused 9 /im behind (A), on (B), and 9 /xm before (C) the track. 
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Sb20 3 , (weight percent), and a refractive index of 1.516 
at k= 0.5467 /xm. A threshold output energy of approx- 
imately 0.5 J was required to produce visible fracture 
damages in this glass. The laser operated in a single 
mode by means of a transmission sapphire etalon 
and had a cryptocyanine-methanol (^-switch. A ballistic 
thermopile-calorimeter was available to monitor the 
integrated energy of each pulse. 




3. Refractive Index Profile 

The "invisible" subfracture filamentary damage 
track in the glass specimen was investigated with a 
micro-Mach-Zehnder twin beam interferometer [3]. The 
glass sample containing the track was ground and pol- 
ished to within optical reach of the microscope objec- 
tive, the reference beam of the interferometer was 
directed through a simultaneously ground and polished 
damage-free specimen of the same thickness. Obser- 
vations were made in the light of the green mercury line 
(Ao = 0.546707 /xm). The correction prisms in the 
reference beam were adjusted until the interference 
fringes were perpendicular to the track. Photographs 
were taken on high-speed Polaroid film (3000 ASA) 2 
and these were rephotographed on line film [4] to 
enhance contrast and eliminate background. Figure 2 
shows the interference photomicrograph of the fringes 
across the track. Some stress birefringence was 
observed when the track was illuminated with polarized 
light. The ray polarized parallel to the track appeared to 
be retarded from the ray polarized perpendicular to 
the track by an amount corresponding to approxi- 
mately 30 percent of the total index change. A cor- 
responding difference in peak displacement was ob- 
served when parallel and perpendicularly polarized 
light was used in the Mach-Zehnder interferometer. 
Since this refractive index anisotropy is a result of the 
cooling history of the track and is difficult to analyze 
accurately, we will consider only the average index 
change in unpolarized light to estimate the minimum 
temperature which the track reached when heated by 
the absorption of laser energy during the self-focusing 
process. 

The fringes show a simple displacement as would 
be produced by a homogeneous cylindrical inclusion 
of refractive index different from the surrounding 
matrix. The shape of the displacement gives no 
evidence of a tubular (hollow) inclusion as has been 
observed by Davit [2]. The base width of the fringe 
displacement measures 3 /xm and it is assumed that 
this is the diameter of the cylindrical track. On its 
apex, the fringe displacement measures 0.15 fringe 
spacings. The direction of the fringe displacement 
indicates that the inclusion has a lower refractive 
index than the matrix. The direction of the displace- 
ment was tested by observing the movement of the 
fringes while blowing a gas of high refractive index 
(freon) through one beam of the microscope. This is in 
contrast to Davit's observation of an increased index 
region. 




FIGURE 2. Photomicrograph of the damage track taken with a 
micro-Mach-Zehnder twin beam interferometer at X= 0.54607 pm. 

The black bar shows 10 /x.m. 



Assuming the track to be a homogeneous cylinder 
of diameter d= 3 /xm and refractive index n 2 embedded 
in a matrix of refractive index n\ , allows calculation of 
the index change in terms of the maximum fractional 
fringe displacement /= 0.15 and the wavelength of the 
light used (in vacuum or air): 



71,-712= An=k fld=2.8 X 10" 



(1) 



2 Certain commercial materials are identified in this paper in order to specify adequately 
the experimental procedure. Such identification does not imply recommendation or endorse- 
ment by the National Bureau of Standards. 



The reduced refractive index indicates that the glass 
has been heated above its transition temperature and 
has been subsequently cooled rapidly. The arresting 
of an elevated fictive temperature structure is well 
known in glasses but under normal quenching condi- 
tions it has not produced such large changes. The 
conditions leading to the reduction in refractive index, 
Ara, are discussed later. 

The lower refractive index, presumably correspond- 
ing to a lower density and to permanent volume expan- 
sion in the cylindrical region, has not caused a sharp 
radial compression in the surrounding matrix. Such 
a sharp compression would be visible as a counter- 
deflection of the fringe. The corresponding compres- 
sion of the surrounding matrix occured apparently 
over a great distance and the associated fringe shift 
is therefore not observable. 
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4. Density of Filament 

The reduced refractive index of the cylinder indi- 
cates that the glass in the cylinder has a density (p 2 ) 
different from that of the matrix (pi). The ratio p 2 /pi 
can be estimated from the Lorentz-Lorenz equation: 



^=(nS + 2)(nS-l)l[(nS 



■l)(/i 2 2 + 2)] (2) 



since n u the refractive index of the undamaged glass, 
is known (n l = 1.516) and A/i has been determined 
interferometrically. p 2 /pi is calculated to be 0.95. 
In other words, the glass in the track has a density 
5 percent below or a specific volume 5 percent above 
the surrounding matrix. A volume expansion of 5 per- 
cent corresponds to a linear expansion of 1.6 percent. 
This figure is used to calculate the heating stresses in 
a later part of this paper. 

5. Thermal History of the Filament 

The density and refractive index of an under- 
cooled liquid such as a silicate glass are dependent on 
its thermal history, i.e., the rate at which the glass is 
cooled from its molten state. The thermal expansion 
coefficient is higher for a liquid than a solid. If a melt 
is cooled rapidly it contracts first as a liquid, but as 
its viscosity increases, the melt is unable to readjust 
its structure in pace with the changing equilibrium 
structure; it solidifies to a glass and then contracts as 
a solid. In normal technological practice, glasses are 
cooled at rates of approximately 2-5 °C/min. At such 
rates the transition from the liquid to the glassy state 
occurs near the equilibrium strain-point temperature 
which is approximately 470 °C in the glass under investi- 
gation. Since the cooling of glasses arrests a high- 
temperature structure, one speaks of the glass as 
having a certain fictive temperature. Figure 3 shows 
idealized curves for a fast and slow cooling treatment 
of a glass sample. 
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Figure 3. Idealized cooling curves of a glass. 

A and B indicate room temperature densities of slowly and quickly cooled glasses 
respectively. 



The ability of a liquid to adjust its structure with 
temperature changes is determined by its volume 
relaxation time. At high temperatures and low viscos- 
ities the relaxation times are short and a very rapid 
cooling rate is required to freeze-in very high tempera- 
ture structures. Consequently, in order to obtain a 
certain fictive-temperature structure, two conditions 
are necessary: (a) the glass must have been heated 
to at least this fictive temperature and (b) it must have 
been cooled at a rate sufficiently fast to prevent re- 
arrangement of the structure when cooling below this 
temperature. 

The index change in the damage track allows us to 
estimate its fictive temperature from the knowledge of 
the thermal expansion coefficients of the glass. We 
also can calculate the instantaneous cooling rates of 
the filament as a function of temperature, compare this 
rate with the relaxation times at these temperatures, 
and estimate the temperature (T/) at which structural 
rearrangements are arrested. Thus, we have two inde- 
pendent means of estimating the fictive temperature. 
The only information to be gained by this, however, is 
the estimate of a minimum temperature (fictive temper- 
ature) which the track reached. Much higher tempera- 
tures may have been reached during the self-focusing 
process, but, for the reasons given above, we only 
have a record of the fictive temperature. 

Details of the calculation of 7/ are shown in Parts 1 
and 2 of the Appendix. We find from the fossil refrac- 
tive index n 2 and the corresponding density ratio, 
p 2 /pi, that the track has at least experienced a temper- 
ature of 1450 to 2800 °C. (The range is caused by un- 
certainty in estimating a value for the high tempera- 
ture expansion coefficient.) In all probability the 
temperature reached was much higher. Since, how- 
ever, the cooling rate of the filament is 10 10 °C/s (see 
Part 2a of Appendix) a fictive temperature not higher 
than 1600 to 1800 °C can be frozen in. At this tempera- 
ture the relaxation time becomes so long (see Part 2b 
of the Appendix) that this high-temperature structure 
becomes arrested. 

An additional clue is provided by the observation of 
microscopic bubbles in the track which indicates a 
minimum estimate of 1500 °C. The manufacture of 
silicate glasses requires that the melt be maintained at 
approximately 1500 °C for many hours in order to 
facilitate the expulsion of dispersed and dissolved 
gases (N 2 , 2 , C0 2 , SO>). At the elevated temperature, 
the melt still contains dissolved residual gases in 
concentrations which correspond to the solution 
equilibrium at the melting temperature and atmos- 
pheric pressure. When such melts are cooled into the 
glassy state, uptake of more gases is minimal because 
of the slow diffusion. In order to produce gas bubbles in 
glass, it is therefore generally necessary to heat it 
above the temperature at which the melt was equili- 
brated with the atmosphere. Under the volume-restrain- 
ing condition of laser track formation, the temperature 
of bubble formation would, however, be considerably 
higher than at atmospheric pressure. 

All the facts discussed above support the assumption 
of considerable heating during track formation — at 
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least to temperatures at which the glass has a viscosity 
below 10 2 P, i.e., the viscosity at which bubble forma- 
tion is expected. This assumption is also supported by 
evidence from other workers. Harper [5] observed the 
extrusion of glass fibers from laser fractures. Hercher 
[1] observed that black body continuum radiation was 
emitted during the damage and estimated a maximum 
temperature of 5000 K. 

6. Thermal Stresses 

Using the model of an expanding infinitely long 
cylinder with a diameter of 3 fjum embedded in an 
infinite solid of the same mechanical properties, we 
have attempted to calculate the stresses in and 
around the damage filament. In the equations, we 
assumed that the cylinder expanded by 5 percent 
of its volume or 1.6 percent of its length isotropically. 
This is the total expansion indicated by our observation 
of the cold, fossil track. The expansion of the hot 
cylinder during the actual laser damage was undoubt- 
edly much higher. We calculate (see Appendix Part 3) 
the minimum stresses in and around the cylinder 
produced during heating of the damage area to be as 
follows: 

A. Outside the cylinder at the surface of the track, 
r=1.5/xm, following the convention that tensile 
stresses are positive, and compressive ones negative, 
the radial, o- rr ; tangential, cree; and longitudinal, 
a Z z, stresses are listed below: 

cr rr = -6.5 X 10 8 N/m 2 = [-95,000 psi] 
o-00 = + 6.9X 10 8 N/m 2 = [+ 100,000 psi] 
<t zz = + 0.46 X 10 8 N/m 2 = [+ 6,800 psi]. 

B. Inside the cylinder, also at the surface, r= 1.5 /xm: 

cr' rr = - 6.5 X10 8 N/m 2 = [-95,000 psi] 
(T 1 ^ =-6.5 X 10 8 N/m 2 = [-95,000 psi] 
(t' zz = - 12.8 X 10 8 N/m 2 = [- 190,000 psi]. 

The tangential stress in the matrix is tensile, and has 
the highest magnitude. Thus, it would preferentially 
lead to a fracture. Its actual numerical value, even 
considering the minimum volume change, is more 
than i of the stress considered to be the inherent, 
internal strength of glass in absence of defects [6]. 
In the presence of the usual defects and considering 
the higher hot stresses around the track, the tangen- 
tial stresses would be sufficient to initiate cracks. 
Such cracks would be expected to occur in planes 
perpendicular to the tangential stress direction, and 
thus would be in half -planes parallel to and containing 
the axis of the track, randomly distributed around the 
circumference. 

The microscopic study of the numerous tracks pro- 
duced in the glass, as described earlier, reveals two 
types of fractures: (a) gross starburst, usually not more 
than one or two per track and (b) very numerous, 



small planar circular fractures along all tracks, associ- 
ated both with bubble regions of the track, and with 
the transparent track region between the bubbles. 

All the planar fractures are almost circular, rest 
with one edge against the track and are in planes con- 
taining the axis of the track. Figure 4 is a sketch of 
the described fractures. We did not find any fractures 
perpendicular to the track as had been described for 
laser damage in polymethyl-methacrylate [7]. We be- 
lieve that the gross starburst cracks were caused by 
foreign inclusions, while the small ubiquitous planar 
fractures are associated with an intrinsic self -focusing 
damage process. 

The morphology of the observed planar fractures is 
in agreement with the findings of the stress calculation 
which predict the dominance of tangential stresses. 

7. Summary 

(1) Laser induced self-focusing damage tracks in 
glass were produced by passive irradiation. 

(2) The refractive index distribution in the damaged 
glass was investigated by double-beam interferometry. 

(3) Filamentary damages were found which consist 
of homogeneous cylindrical regions, 3 /xm in diameter, 
having a refractive index of Arc = 2.8 X 10~ 2 below the 
surrounding glass. 

(4) The reduced refractive index in the damaged 
region is calculated to correspond to a reduction in 
density by a factor of p 2 /pi — 0.95. 

(5) It is calculated that the glass in the filament has a 
structure corresponding to a Active temperature of 
1450 to 2800 °C, the range being caused by uncertainty 
in the value of the high temperature expansion 
coefficient. 

(6) It is calculated from heat flow conditions and a 
typical relaxation spectrum of glass that it is not pos- 
sible to freeze-in a high-temperature structure corres- 
ponding to a higher Active temperature than 1600 to 
1800 °C, even though the actual temperature during 
laser beam exposure may have considerably exceeded 
this temperature. 




Figure 4. Sketch of the small planar fractures observed along 
all tracks. 

The fractures are in half-planes containing the track axis (z) and a random radial 
vector R. 
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(7) Bubble formation in damage tracks is evidence 
that the laser beam heated the glass at least to 1500 °C. 

(8) Equations for the axial, radial and tangential 
stresses in and around a heated cylindrical region in 
glass are derived. It is shown that for heating the cylin- 
der to its Active temperature only, corresponding to 
a volume change of 5 percent, the tensile tangential 
stress reaches as much as 100,000 psi which is i of 
the inherent, internal strength of glass in absence of 
defects. 

(9) A microscopic investigation o£ the laser tracks 
shows large numbers of small planar fractures originat- 
ing from one side of the track and oriented parallel to 
the track. Such fractures have the proper orientation to 
be caused by the predicted dominant tangential 
stresses arising from laser-beam heating. 

(10) It is concluded that thermal expansion of the 
glass is a sufficient mechanism to produce fractures 
within the glass. This conclusion is supported by the 
fact that glasses with low thermal expansion (fused 
silica, borosilicate glasses) show a higher damage 
threshold to passive laser beam exposure than glasses 
with high thermal expansion [8]. 

8. Appendix 

Parts I and II of this appendix show detailed cal- 
culations of the estimated ranges within which the 
Active temperature of the damage track occurs. 
These calculations rest upon assumptions and equa- 
tions which are well followed, in general, by silicate 
glasses, and are only proposed to give an estimate of 
the minimum temperature reached by the damage 
track. Each part presents totally independent analyses 
which rest on different assumptions and available 
data. 

Part III is a rigorous analysis of the stresses involved 
in the heating process and does not depend upon the 
results of Parts I and II. Since we are calculating stres- 
ses developed in heating, and the heating process was 
very fast, we may neglect the effect of plastic defor- 
mations, or volume relaxations and compaction effects. 
Therefore, Part III is only concerned with the elastic 
problem of heating a cylindrical filament within a 
cold matrix. 

Part I. Calculation of Fictive Temperature from 
Index Change 
The estimate of a fictive temperature from a change 
in index depends upon the variation of index with 
specific volume, and of specific volume with tempera- 
ture. From the Lorentz-Lorenz equation, we can easily 
estimate Av and An, but the dependence of specific 
volume on temperature up to high temperatures is 
only known in general terms. We know the expansion 
coefficient of this glass up to 300 °C and may apply it 
as high as the transformation range (550 °C). We also 
have measured the temperature dependence of index 
in the transformation region from 550 to 580 °C. 
From these data, we may calculate an upper and a 
lower temperature estimate for 7/, as shown in figure 
5. The relative size is plotted versus temperature 
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FIGURE 5. Extrapolation procedure for estimating the high 
temperature specific volume of the glass. 

T\ and Tt are the upper and lower estimates of the fictive temperature. The measured 
solid expansion coefficient is shown. The expansion coefficient quoted for the transforma- 
tion range (550-580 °C) has been calculated from the index data by using the Lorentz- 
Lorenz equation, and is only shown here for illustration. 



and follows the general behavior expected of silicate 
glasses. The estimate, T\, is an extrapolation of the 
index dependence on temperature. It uses the expan- 
sion coefficient characteristic of the transformation 
region, and will underestimate Tj. The estimate T>> 
uses the solid expansion coefficient and will overesti- 
mate Tf. 

Measurements of index change as a function of 
temperature in the transformation range from 550 
to 580 °C are shown below: 



Table 


1 


/i(A = 0.5895Vm) 


Temperature 


1.51337 
1.51363 
1.51406 
1.51436 


584 °C 
575 °C 
562 °C 

552 °C 



These points plot as a straight line whose slope yields 
the temperature dependence of index: 



W\- 



3.16xiO- 5 /°C 



The observed change of An = —0.028 gives a tempera- 
ture change of 900 °C. The base glass is assumed to 
have a fictive temperature of 550 °C, therefore we 
estimate 7\ - 1450 °C. 

In the calculation of 7^, we now invoke the Lorentz- 
Lorenz equation to calculate the difference in volume 
between the track and matrix. If n\ is the index of the 
matrix at 1.516, corresponding to a matrix fictive 
temperature of 550 °C, and An is — 0.028 we may cal- 
culate n 2 to be 1.488 and calculate the size ratio 
L 2 /Li to be 1.016. The solid expansion coefficient is 
7.15 X 10 _6 /°C, therefore, the corresponding tempera- 
ture change is calculated to be 2240 °C. Adding the 
fictive temperature of the matrix leads to an estimate 
for Ti — 2800 °C. The fictive temperature estimate is 
therefore seen to be between 1450 and 2800 °C by 
this method. 
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Part II. A. Calculation of Cooling Rate of Filament 
The hot filament may cool by thermal conduction 
and/or radiation. Due to the good thermal contact 
assumed between the hot zone and its surroundings, 
the heat flow by conduction is much larger than by 
radiation, and the latter may be totally neglected. 

The problem of a hot cylinder embedded in a matrix 
has been solved by Carslaw and Jaeger [9] for the 
change in temperature of the matrix as a function of 
time and distance from the hot zone. Due to the good 
thermal contact between the cylinder and matrix, we 
expect a very high cooling rate, therefore, we are only 
interested in the short-time solutions. The cooling 
rate may be calculated from the heat flow between the 
hot cylinder and a larger concentric cold sink surface. 
We shall, therefore, use the solutions of Carslaw and 
Jaeger to estimate the distance of this cold sink from 
the hot zone. An analysis of their calculations [fig. 41, 
p. 337, ref. 9] leads to table 2 below, assuming a heat 
conductivity, K= 0.0121 J cm -1 s _1 °C _1 ; a specific 
heat, C= 1.34 J g _1 °C _1 ; a density, p — 2.1 g/cm 3 , and a 
filament radius a= 1.5 (xm. 

Table 2 







Time taken to increase 


Distance from filament 


Temperature at / () : 


the temperature of r () 


ro/a 


7yr(filament) 


to the value To 


10 


0.4 


8X 10-*s 


10 


.2 


8x 10- 3 


1.4 


.4 


3.2 X10- 5 


1.17 


.4 


2.4X10- 6 


1.17 


.2 


8X10- 7 


1.09 extrapolate 


.4 


8x 10- 7 


1.04 extrapolate 


.4 


2.4 X10- 7 



Since we are interested in finding the cooling rate 
applicable for the first 200 °C drop in temperature, and 
the cooling rates we expect to obtain are of the order 
10 9 to 10 10 °C/s, let us consider a cold sink at r = 1.17a 
for a slow rate, and r = 1.04a for a fast cooling rate 
estimate. The cooling rate can be calculated as 
follows: 

dT = 1 dq = 27raLK(T m -To) = 2K(T n -To) 
dt pVC dt C irpa 2 L(ro — a) C pa (ro — a) 

where K is the thermal conductivity, dq/dt is the rate of 
heat flow, C is the specific heat, V is the filament 
volume, L is a unit length and To is the matrix tempera- 
ture at r . The low and high rates (r = 1.17a and 
r = 1.04a) are calculated respectively to be: 

(^)=2.3xionr fi ,-r ](°c/s) 
(^)=9.5xionr fil -r ](°c/s). 

The resulting heat conduction cooling rates are listed 



in table 3 below and it is seen that the filament cools 
fast enough to undergo a temperature change of 
approximately 100 °C in less than 4 X 10 -8 s, and there- 
fore our assumed cold sinks remain so during that 
period. In the calculation of cooling rates, it is, there- 
fore, suitable to assume that To represents room tem- 
perature, 20 °C, and thus ignore it. 





Table 3 




Filament temperature 


\ dt )„ 


\ dt ), 


1000 °C 


i.oxio io x:/s 


2.4X10 9O C/s 


1200 °C 


1.2 




2.8 


1400 °C 


1.4 




3.3 


1600 °C 


1.6 




3.8 


1800 °C 


1.8 




4.2 


2000 °C 


2.0 




4.7 


2200 °C 


2.2 




5.2 



Part II. B. Calculation of Volume Relaxation Time 
for the Filament 

The volume relaxation time at high temperatures 
may be estimated by comparing this glass to a sodium- 
borosilicate for which high-temperature data have 
been reported [10]. At temperatures above 1400 °C 
most oxide glasses are found to have a single or narrow 
distribution of relaxation times, therefore, character- 
izing the most probable relaxation time will suffice to 
analyze the relaxation processes. Results from the soda- 
borosilicate study [10] also show that at high tempera- 
tures the volume relaxation time is approximately 3 
times the shear relaxation time. The latter may then 
be found by considering the shear viscosity data, and 
comparing it with the soda-borosilicate (SBS) glass [10]: 



-=-2L 

T?SBS 



T SBS- 



This equation rests upon the likely assumption that the 
moduli of the two glasses are equal. This falls well with- 
in the order of magnitude estimate which we are seek- 
ing. The volume relaxation time is thus found as: 



r v = S 



_2_ 

^SBS 



' T SBS- 



Once r v is calculated, we have the rate at which the 
structure rearranges itself. The minimum cooling rate 
necessary for freezing a given structure is found by re- 
quiring t v to be the time necessary for the relaxation 
time to increase by a factor of 2. This factor approxi- 
mates the width of what is usually assumed to be a 
single relaxation process. Thus, by the time r v has 
elapsed, it would require the structure 2 r v to rearrange 
and so on until room temperature is reached. Once the 
cooling rate necessary to freeze a structure equals the 
heat conduction cooling rate, then the latter is a likely 
mechanism for arresting structural changes at the tem- 
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perature of coincidence. The volume relaxation time 
and minimum cooling rates are given in table 4 below: 

Table 4 



Temperature 


M 


nimum cooling rate 


\ olume relaxation 










time T r 


1000 




2.5X10 <;C 


C/inin 


10 5 s 


1200 




10 8 




5X10- 7 


1400 




10 9 




3 x 10- 8 


1600 




8X10 9 




6 X 10- 9 


1800 




4X10 10 




2 x 10- 9 


2000 




1.5 X10 11 




5X10- 10 


2200 




4X10 11 




2X10- 10 



Referring to table 3 we may now compare the two 
aforementioned cooling rates, and find the point of co- 
incidence. An upper and a lower estimate of the con- 
duction cooling rate were calculated, therefore, we have 
two temperature estimates. The temperatures ob- 
tained from these calculations give us an estimate of 
fictive temperature from 1600 to 1800 °C. 

Part III. Solution of Stresses Developed by Heating 
of the Laser Track 

In this section, we are interested in calculating the 
stresses developed both inside and outside the cylinder 
representing the laser track when its temperature 
reaches a value corresponding to the measured change 
in index. By the Lorentz-Lorenz equation, we can 
calculate the change in volume associated with the 
index change without having to speculate about the 
temperature dependence of the index of this glass. 
The stresses can then be directly calculated from the 
volume change. 

After a somewhat detailed look through the litera- 
ture, for a general solution of the problem, we un- 
fortunately found that the most common solutions 
involve the "plane-strain" approximation. This ap- 
proach assumes that there is no strain in the longi- 
tudinal direction and that the necessary stress is 
provided externally to the system. The solution does 
not represent the problem at hand, since a large 
longitudinal stress is expected to develop in the cold 
matrix. The need to determine all stresses in the matrix 
(outside the hot cylinder) requires a more general 
solution which is given below. 

Some solutions of the generalized plane-strain 
problem may be found in the literature [13], where, 
in this case, the longitudinal stress is assumed to be 
a constant. We have not found one which proceeded 
beyond the first motions of solving the problem, 
therefore, we have chosen to solve the problem from 
the general stress-strain relations. 

The problem is best expressed in cylindrical co- 
ordinates, since the angular symmetry eliminates the 
^-dependence of the solutions. The general stress- 
strain relations used here are taken from Love [11], 
and particularly from Boley and Weiner [12]. 

The solutions derived below are based upon the 
assumption that the positive volume change of the 
cylinder in the longitudinal (z) direction occuring during 



the heating process, induces stresses outside the 
cylinder which occur throughout the cold matrix and 
reach a maximum at the interface. Outside the 
cylinder, the longitudinal extension (U z ) is therefore 
proportional to the length [or distance away from the 
origin which we locate in the middle of the cylinder], 
and the longitudinal strain is uniform. The latter, 
€**, is only a function of the radial distance: 

U z =f(r)z; tzz=f(r) 

This allows solution of the equilibrium equation in the 
axial direction: 





d&rz ()(Tzz , (Trz _ ~ 

dr dz r 


where 






CTrz = /JL 


~dU r dUj 

dz dr 


and 







(Tzz= (a + 2/x) € zz +k€rr+teee 
a and /jl are the Lame elastic constants for the material, 



and the strains, c u - \ 
dij are defined in Love 



and the stresses 



m_ dUj 

dXj dXi\ 

11 J, and Boley and Weiner 
[12]. The cylinder is considered infinite to avoid end 
effects, therefore, the radial displacement, f/ r , has 
no 2-dependence. Since neither e 2Z , nor e rr , nor €00 
depend on the z-variable, the differential equation 
is easily solvable for/(r): 

cPf(r) | 1 df(r) =Q 
dr 2 r dr 

and leads to the following solution with Bo and ro as 
integration constants. 

/(r)=fioln(r /r) 

Uz = So 2 In (r /r) 

e Z z = Bo In (r /r). 

The equilibrium equation in the radial direction may 
now be written for the region outside the cylinder: 



dr 



■ + 



d(Trz , (T n 



dz 



(TOO 



= 



ith: 



(Trr — ( A + 2/Ji) — hX h A — — 

dr r dz 

(Tee = (A + 2/x) h a — h a — — 

r dr dz 



(Trz — — fJiBoz/r. 



The resulting differential equation is again easily 
solvable: 
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im^Hm^ 



with the following solution: 



U r (r)=- 



Bo 



X + 



X + 2/a 



g]"" 



(ro/r) + 0/r. 



The characteristic length of the logarithmic term 
ro, is taken to be the same as for U g . While this assump- 
tion is made for simplicity in solving the problem, 
it is reasonable, as ro characterizes the penetration of 
the effect of the longitudinal stress into the matrix. 
Q is a constant whose value will be determined by 
the boundary conditions. The logarithmic term in 
both solutions is only valid up to r=ro. The following 
expressions are then obtained for the stresses outside 
the cylinder: 

CTrr — (A + fl) Bo/2 



(Tee 



-m&] 



Bo/2 



[-&&>• 



In (ro/r) + 2/u <?/r 2 



Inside the cylinder, the solutions resemble the 
"plane-strain" approximation with the exception of 
an allowed expansion of the cylinder in the longitudinal 
direction [see ref. [12], pp. 288-291]. 



U'=Az 



3A. + 2fx 



2ix_ 



m 



AL 



is the total fractional expansion of the material 

in heating to the temperature of the track, and A is 
a constant to be determined by the boundary condi- 
tions. Since the cylinder is small, we have not allowed 
a variation of U z with radial distance. This approxi- 
mation proposes that the rod expands uniformly in 
the longitudinal direction, and greatly simplifies the 
problem without much effect on the final conclusions. 
The resulting stresses may be written as follows: 



cr r 



, = _ J >(3X + 2/i)1 AL 



L (A + 2//,) J L 



, _ r /x(3X + 2/x) l AL 

o*-™ [ (X + 2/t) J L 



<r' zz = (\ + 2h)A 



= A. 



r 2/x(3A + 2 i ic) 
L (X+2/i) J 



AL 
L 



We now turn to boundary conditions to solve for 
the four unknown constants: A, Bo, ro, and Q, assuming 
the following typical values for the Lame constants 

and the calculated elongation, -=- : 



A =2.1 Xl0 10 N/m 2 

/* =2.66X10 10 N/m 2 

AL/L = 0.016. 

The boundary conditions are as follows: 

(i) the radial stress is continuous across the 
boundary: 



CTrr (a) = Or' r r(a) 



B.C.-l 



(ii) the tangential strain and the radial displacement 
are continuous across the boundary: 



U r (a) = U' r (a) 



B.C.-2 



(iii) the longitudinal strain is continuous across the 
boundary: 



e**(a)= €«(a) 



B.C.-3 



(iv) there are no external longitudinal forces on the 
sample: 



ds = 



I (T zz 'ds+\ 

J So J Si 

where So is the cross-sectional area. This becomes: 
a zz rdr+ cr' zz rdr = 

Ja JO 



B.C.-4 



Starting with B.C. -2 and B.C.-3 we calculate Q 
from an estimate of A : 

<? = 2.82X10- 14 + 7.23X10" 13 ^. 

Then, Bo and ro are calculated in turn from B.C.— 1 
and B.C. -3 as follows: 

£ = 1.29A + 9.95 XlO n Q- 2.8 XIO" 2 

r =EXP^- 13 - 41B °" 






Bo 



']■ 



Finally, a new value of A. is computed from B.C.-4 
and reinserted into B.C.-2. The integrations of B.C.-4 
for the region outside the cylinder are limited to ro 
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since the solutions are only valid in this range. They 
lead to: 

A = lO" 2 + 2.78 X 10 - 5 flo ~ 1.23 X 10 12 fl o rg. 

The iteration converges quickly and the following 
values are obtained for the constants, after we reached 
a variation smaller than 0.5 percent of each term: 

A =6.22X10" 4 

Q =2.86XlO" 14 m 2 

Bo=1.31X10- 3 

r =2.46X10" 6 m. 

The stresses at the surface are then computed for 
r=a= 1.5 fim. 

OUTSIDE" 

o-rr =-6.5 X 10 8 N/m 2 = - 95,000 psi 

o-ee =+6.9 X 10 8 N/m 2 = + 100,000 psi 

o-^= + 4.6 Xl0 7 N/m 2 = + 6,800 psi 

INSIDE 

o-' 2Z = -6.5 Xl0 8 N/m 2 = - 95,000psi 

cr'^ = -6.5 X10 8 N/m 8 = - 95,000psi 

o-' zz = - 12.8 X 10 8 N/m 2 = - 190,000 psi 

Tensile stresses are positive while compressive ones 
are negative. Inside, the radial and tangential stresses 
are equal as in the plane strain problem, but their 
magnitude is slightly reduced due to the allowed 
longitudinal strain through the shear modulus. The 
longitudinal stress is seen to be very large, and all 3 



stresses are compressive, as expected. Outside, the 
tangential and longitudinal stresses are tensile, with 
the former equaling three times the latter. The low 
value of the external longitudinal stress at the interface 
stems from the elastic distribution of the total force 
far into the matrix. 
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